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^ 1 Introduction 

O 

The problems in plasma physics always require some analysis of the collisions. An equation 
which has a general form describing collisions is the "Boltzmann Equation". It has the form 
as follows: 

Where f a (x,v,t) is the Boltzmann function of a particle in the usual notation, a is the 
macroscopic contribution to the acceleration of the particle at (x, v, t) and the term on the 
right represents collisions. A commonly used form of the collision term is of the Fokker- 
ft Planck type[l] 

Here A^ and B^ are the "friction" and "diffusion" coefficients respectively. Several different 
kinds of forms of At and Bij are provided by previous plasma physicists. A widely used 
kind of collision term is called "Rosenbluth Collision Term" [2], which can be written in the 
• <-H following form: 

ft [ dt )c ~ 4- al3 dv' [Ja <9v 2dv' dvdv } [ ) 
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In this collision term, H a p and G a p are called Rosenbluth potentials. They satisfy following 
relations: 

ON Ha/?(v) = / dv'-: (4a) 

mp J |v-v'| 

GW(v) - | dv'//j(v')|v-v'| (4b) 
2 V^ = - 47r(m " +W ^ /. (4c) 

Several other researchers have paid much attention about the eigenvalues and eigenfunctions 
of this collision term. Lenard and Berntein(LB)[3] reduced the collision diffusion coefficient 
to be a constant related only to thermal velocity [Eq. |(5|]. Then they got an exact analytic 
solution with a dispersion relation. After their work, C.S.Ng,A.Bhattacharjee and F. Skiff 4 
further discussed the completeness of LB collision term theoretically. Their results indicated 
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that the LB operator has a discrete eigenmode spectrum, totally different from Case- Van 
Kampen eigenmodes[5] [6]. 

i-W)o = ^Ma + v ^) (5) 

J. P. Dougherty [T] derived his own kinetic equation based on the conservation theorems [Eq.Q]. 
Then M.W.Anderson and T.M.O'Neilflj analytically gave a theoretical formulation of all the 
eigenvalues and eigenfunctions of the Dougherty collision term. 

/ 9 fa \ 9 KTdf a 

However, the previous researchers could only deal with Fokker-Planck collision term by re- 
ducing its complexity and difficulty to a much simpler mathematical form. For a general 
case, there are still no thorough and detailed discussions. 

In this paper, we utilize Rosenbluth collision term for only ion-ion collisions and discuss 
the eigenvalues and eigenfunctions of this operator. Section 2 mainly focuses on the math- 
ematical deducing procedure of linearizing Rosenbluth collision term. By comparing the 
whole operator and the differential term in the operator, Section 3 discusses the minimum 
eigenvalues of the least damped wave modes and their corresponding eigenvalues. Section 
4 concentrates mainly on the theoretical analysis of the differential term including a proof 
of the completeness of the eigenvalues and orthogonality of the eigenfunctions calculated 
only by this term. Section 5 gives further discussions about the differential term. Some 
mathematical transformations are used and the continuum of the corresponding eigenvalues 
spectrum is indicated. Section 6 sums up the all conclusions. 



2 Linearization of Rosenbluth Collision Term 

For research convenience and experiment reality, Fokker-Planck equation with this collision 
operator is always linearized to its first order perturbation by assuming the zeroth distribu- 
tion function as the a Maxwellian form / M [5]. 

/ = f M + f l e< kz ~ ut) (7a) 
E = -A<p, where <j> = <j)o + fae i{ - k *- ut '> (7b) 

After that procedure, one can get: 

(^) c = i(kv z - co)fl + % Z -^pLfM = c [f M, f] + C[f\ n (8) 

The linearized collision operator includes two parts: C[f M , J 1 ], integral part and Cf/ 1 , / ], 
the differential part. These two parts can be written in the form: 

c[/\n = -^-(/'^-^-/^) W 



The collision terms in Eq.(9a) and Eq.( |9b[ )can be simplified by working in spherical coordinates:^ = 
|v|,/x = v z /v,and by using symmetry properties of f M and f 1 , 

Coif =C[f J]=r[^^) (4tt/ / -— / ) (10) 
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1 (11) 

(H i, Goi) an d (Hi , G\o) refer to Rosenbluth potentials of perturbation f 1 and background 
f M , respectively. The collision strength parameter T,Hiq and Gio are given by 

_ 4vrZ 4 e 4 ; 3^ 4 

T = ^ — /nA = — ^ (12a) 

m| An 

_ 2(27r)3/a ^e-" a ^ -er/(,/^) 



<9v 

(27r) 3/2 er/(i;/\/2) (12c) 
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<9?; 2 dv 

Where erf(x) is the error function. To make analysis more convenient, we simply make the 
normalizations, f M — > f M [ni/(2ir) 3 / 2 vf] 2 and f 1 -» f 1 [n i /(27r) 3 / 2 i' 3 ] 2 and denote C w as the 
differential part of the operator and Goi as the integral part. 

in reference [S]and[S], the perturbed ion distribution function f 1 is written as a sum of 
Legendre polynomials Pi, 

f 1 = "£ai(v)Pi(v) (13) 

1=0 

The expansion by using Legendre polynomials is also called "Moment Expansion" .If there is 
no electric field and the perturbation is isotropic both in position space and velocity space, 
the equation of individual a; will be decoupled. This makes analysis of collision term much 
more easier, and that is why moment expansion is so widely used. 



3 Eigenvalues and Eigenfunctions 

For further discussion, here we only consider f 1 to be isotropic in position place and velocity 
space and all the calculations are under the assumption of no electronic field. Then the 
Fokker-Planck equation becomes: 



-3y/TT d v 8H W 



v 2 k\ dv 2 dv 
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(14) 



And Hi and Gj are Rosenbluth coefficients: 



H,(v) 



d^jv) 
dv 2 



-47T 
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' ai (v')v' 2+l ( l{l - 1] {l + m + 2){l - 1/2)U ' 2 ' 



(l + 3/2)v l + 3 



dv' 



+ 00 



ai(v') 



„H-3 



(Z + 3/2)v' : 



(15a) 
(15b) 



Eq.(14)is the eigenvalue equation. The first term in { } in Eq.(14) is form the differential 
part and the second term in [ ] comes from the integral part. 

In order to calculate the eigcnfunction A of Eq. ( 14 ) , a common method by expanding ai in 
Sonine polynomials S„ +1 (u 2 /2) has been discussed by previous researchers [8] [9] . But in our 
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work, numerical simulation is employed by finite diffcrncing the derivative terms and dis- 
cretizing the integrals on velocity grid Vi = {v\, V2, • ■ ■ , Vj max }, where Vj max is the maximum 
velocity used in simulation and it is normalized to thermal velocity, jmax is the maximum 
number of velocity girds. By comparison, the eigenvalues of the differential operator are also 
calculated. As to our results, only the minimum eigenvalues are taken into consideration 
because their corresponding eigenmodes are of least resistance. And it is very obvious from 



Eq.(14) that all eigenvalues u> should be imaginary numbers. Therefore, in our discussion, 



only the imaginary parts of eigenvalues are taken into account. The results are shown as 
follows: 
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(a) Figure 1 



(b) Figure 2 
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(c) Figure 3 



Figure 1 plots the case of / = 0,v m = 7vi, the horizontal axis means the number of grids 
in each case. The dots in black represent the eigenvalues calculated by the whole operator, 
while the read ones represent the eigenvalues calculated by merely differential operator. By 
examining figure 1 to figure 3, the eigenvalues of the whole collision term and the differen- 
tial term are almost the same. Especially for large I and v m , their differences become much 
smaller. The oscillations of eigenvalues of the whole collision term may be a result of the in- 
fluence of integral terms. The result that the eigenvalues are approximately the same applies 
to other cases as well. Then we go deeper to examine the eigenfunctions corresponding to 
the minimum eigenvalues calculated respectively by the whole term and only the differential 
term. 
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(d) Figure 4 



(e) Figure 5 




(f) Figure 6 



In figure 6,7,8, the maximum velocity is 8Vi, and the number of girds is 500. The only 
difference is the subscript I. The line with blue color is the corresponding eigenfunction of 
the minimum eigenvalue of the whole term and the green line is of the differential term. 
These two kinds of eigenfunctions differ in small velocity but become more alike in the large 
velocity region. All the eigenfunctions should converge to when v — > +oo by considering 
physical facts. By comparing these 3 graphs, only when I = 0, the eigenfunctions do not 
converge to zero when v approaches 0. For higher order, say Z > 0, all the eigenfunctions 
start from value of zero. The reason of this phenomenon will be explained in following 
sections. 



4 Eigenvalues of Differential Operator 

In this section, we only talk about theoretical properties of differential term in the whole 
operator. The differential operator have been shown in Eq.(ll). After neglecting integral 
operator, the Fokker-Planck equation is like: 



- iojf 1 = Cio/ 1 



To write Eq.(16) in a more detailed way, we get: 

-[(1 



l (2ir) 3 / 2 vf 
1 dG w d 



2v 3 dv <9/i 



_L__<9_ dH 10 _ v 2 
2v 2 dv dv 



(16) 



(17) 
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By writing — iu)[(2w) 3 ' 2 vf /rii] 2 /T simply as A and then using moment expansion[Eq.(13l], 
we gain the eigenvalue equation of the differential operator. 

l{l + l)dG w 1 d r dH w _ v 2 /2 d 2 /2 , 



Eq.(17) describes the evolution of distribution function according to time perturbed by a 
test particle. However, the background remains Maxweillian distribution not affected by the 
test particle. Only considering differential operator costs the equation to lose some basic 
properties of the whole Rosenbluth collision term such as the conservation of momentum 
and energy, even though the number of particles still remains conserved. This problem may 
brings out some unphysical results but in some cases when the differences are not so obvious, 
the results are still reliable and valid. 



Eq.(18) can easily be transformed into Strum-Liouville(SL) type differential equation: 

^ [k(v) ■ X'(v)] - q(v)X(v) + Xp(v)X(v) = (19) 



In Eq.(19l, the specific functions are: 

X(v) = ~e v2 / 2 ai{v) (20a) 

k(v) = -ve- v2 ^^ (20b) 

q (v)= l V±Ve-" 2 ^ (20c) 
V dv 

p(v) = 2v 2 e~ v2/2 (20d) 



The boundary conditions arc: 



lim X(v) < +oo 
lim X(v)=0 

v— >+oo 



(21) 



According Strum-Liouville theorem, when I = 0, it is easy to prove that all these functions 
are positive. Thus, all the eigenvalues comprises a complete set and all the eigenfunctions 
are orthogonal. For orthogonality, the weighing function is p(v). Thus, the completeness of 
the eigenvalues and the orthogonality of the eigenfunctions when I = are proved rigorously. 



5 Further Discussions on Differential Operator 

In order to get a deeper insight into differential operator, first we reduce the complex differ- 
ential equation Eq.(17l to a simpler form like this: 

d 2 ai 



Cwai 

A,B,C are functions of v, 

A = 
B = 
C = 



A' 



dv 2 



ov 



Xai 



1 dH 10 
2v dv 

1 ,dH 1Q 
2v 2 dv 



d 2 H u 
ov z 



,dHu 
dv 



ldH w ld 2 H w 1(1 + 1) dG 
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v dv 2 dv 2 



2v 3 



dv 



(22) 

(23a) 
(23b) 
(23c) 



G 



In order to eliminate the first derivative order term in Eq.(22), namely dai/dv, a mathe- 
matical transformation called "Liouville-Green Transformation" is used: 



a t (v) = H(v) - exp(-~ I L(r'ulr') 



(24) 



Here we donote exp(— | J Q ^(v')dv') as a new operator O, Eq.(24| can be written in another 
way: 

ai — H ■ O = O ■ H (25) 



After some mathematical procedure, the Eq.([22|) can be expressed in the following form: 

(26) 



(O-'C^H = A{H" + [r)(v) - \e(v) - ~?(v)}} = XH 



Where £(v) = B/A,r](v) = C/A. Then, a new operator can be defined: P = O 1 CioO, its 
eigenvalue equations is: 



PH = A{H" + [„(«) - \e{v) ~ \i'{v)]} = XH 



(27) 



By comparing the form of Eq. ( 22 ) and Eq. (271, the eigenvalues are the same for operator do 
and the operator P. And this is somewhat like an orthogonal transformation by carefully 
examining the components of P. 
We also define U(v) as a potential like function. 



U(v)=v(v)-y\v)-~Z'(v) 



(28) 



This potential like function will determine the property of continuum of the eigenvalues. It 
will be discussed later. More over, Eq.(27) are much like a wave equation, and this makes 
the physical essence more clearer. 



5.1 Behaviors of U(v) when v — > 

Because both in £(v) and 77(f), when v — > there may be a singularity, it is very necessary to 
have a thorough discussion about the behaviors of the related functions when v approaches 
0. 

First of all, error function can be expanded using Taylors formula at v = 0: 



2 11 

erf(x) = ^(x--x 3 + —x 5 + o(x 7 )) 
y/n S 10 

Then function dHio/dv can be written in progressing form: 

dHiQ 8 5 o , c. 

— - — ~ irv irv + o[v ) 

ov 3 4 



By using Eq.(12c), we get progressing form of dGio/dv: 

dG 1Q 



dv 



8 16 3 / 5n 
-7T1> irv + O (V ) 

3 15 y ' 



(29) 



(30) 



(31) 



By using Eq.(23a)-Eq.(23c), A, B,C have the form: 



A ~ --7T + 0(V 2 ) 

3 

8 1 

B~ --„(-) + o(v) 

6 V 

, . 81(1 + 1) , 41(1 + 1) ,1, , 2 s 
C ~ -4tt + V 7 7T + \ > n(-) + o(v 2 ) 
15 3 



(32a) 
(32b) 
(32c) 
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Finally, we get the progressing behavior of U(v): 



tf(«)~(3- 



21(1 + 1) 



)-l(l + l)(-)+o(v) 



(33) 



Eq.(33) gives some explanations stated about eigenf unctions in section 3. When it is the case 
that I = 0, the term with a singularity form 1/v 2 disappears. Thus, the eigenfunction should 
not be zero in the v — > limit. But when I > 0, this singularity remains. The potential like 
function is infinite when v — > 0. That is why for I = 1, 2 the eigenfunctions start from to 
avoid infinity. 

This phenomenon can also be viewed in a physical way. When v — > 0,the ion collision 
rate becomes infinite, so the perturbation tends to become more isotropic. Hence, the 
contribution of Pq, which describes the isotropic part of distribution function, becomes the 
largest and the other Legendre Polynomials become less important. So only for I = the 
eigenfunction has a totally different behavior in zero velocity limit. 



5.2 Continuum of Eigenvalue Spectrum 

The continuum of eigenvalues of differential operator can be proved by examining the poten- 
tial like function U(v). Just like potential barriers and potential walls in quantum mechanics, 



if the following relation) 34 1 is satisfied, which means there is no potential barrier to trap 



the eigenmodes, the eigenvalues will be continuous. 

U(v) < (34) 
Here are the graphs for I = 1,2: In Figure 7 and Figure 8, the potential like function 
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(g) Figure 7 



(h) Figure I 



U(v) is negative which indicates relation) 34) is satisfied. The continuum of the eigenvalues 
of differential operator also give some evidence that the eigenvalue spectrum of the whole 
Roscnbluth operator may be continuous as well. 



6 Conclusion 

In this paper, firstly we discuss the general collision term of Fokker-Planck equation. Then, 
by using moment expansion, we get the eigenvalue equation of Rosenbluth collision operator. 
Some numerical results have been presented, along with the minimum eigenvalues calculated 
by only using differential operator of the whole operator. These results point out that in 
some region these two eigenvalue spectrums are almost the same but as to eigenfunctions, 
there are some differences in low velocity region. 
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A new method on analyzing differential operator is proposed. The eigenvalue equation can 
be transformed into Strum-Liouvillc type differential equation. Then the completeness of 
the eigenvalues and the orthogonality of eigcnfunctions when I = are therefore proved by 
referring to Strum-Liouville theorem. 

To further analyze the property of differential operator, Liouville-Grccn transformation is 
introduced to reduce the complexity of differential operator. A new operator P = 0~ 1 CioO 
which has the same eigenvalues is defined and all the problems are from the potential like 
function. Then, we thoroughly examine the behavior of U(v) and get some valid results. The 
continuum of eigenvalues of the differential term is indicated by the fact that the potential 
is always below zero. 

However, some properties of the whole operator are still difficult to interpret because of the 
integral term which is to complicated to analyze theoretically. More over, some conservation 
relations are neglected if only differential term is taken into consideration. Hereby, there are 
still much analytical and numerical work remaining to be done. 
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